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; ABSTRACT 

■ We discuss vortex-mediated mutual friction in the two-fluid model for superfluid neu- 
tron star cores. Our discussion is based on the general formalism developed by Carter 
and collaborators, which makes due distinction between transport velocity and mo- 

| ' mentum for each fluid. This is essential for an implementation of the so-called en- 

trainment effect, whereby the flow of one fluid imparts momentum in the other and 
£■ — ' vice versa. The mutual friction follows by balancing the Magnus force that acts on 

lO . the quantised neutron vortices with a resistive force due to the scattering of elec- 

' trons off of the magnetic field with which each vortex core is endowed. We derive the 

form of the macroscopic mutual friction force which is relevant for a model based on 
1^ ' smooth-averaging over a collection of vortices. We discuss the coefficients that enter 

f^i . the expression for this force, and the timescale on which the two interpenetrating fluids 

' in a neutron star core are coupled. This discussion confirms that our new formulation 

accords well with previous work in this area. 

6 ■ 

■ 1 INTRODUCTION 

C3 ■ 

| A superfluid rotates by forming a dense array of quantised vortices. In the case of mature neutron stars, which are expected to 
• ^ ■ exhibit large scale superfluidity since their core temperatures are orders of magnitude below the Fermi temperatures for both 
r> ' neutrons and protons, these vortices should play a key role in determining the rotational dynamics. In fact, the sudden spin-up 

■ associated with the observed radio pulsar glitches (Lyne, Shemar & Graham Smith 2000), and the relaxation that follows, 
is commonly viewed as strong evidence of transfer of angular momentum between a superfluid component and the charged 
component to which the star's magnetic field is locked and which, presumably, is linked to the pulsar emission mechanism. 
Hence, the notion that a neutron star acts as a kilometer-sized superfluid system is strongly supported, both theoretically 
and observationally. 

The purpose of this study is to model vortex-mediated dissipation in the two-fluid paradigm for superfluid neutron stars. 
The common view is that the most important dissipation mechanism in a superfluid neutron star core originates from electrons 
scattering off of the magnetic fields associated with the individual vortex cores (Alpar, Langer & Sauls 1984; Mendell 1991b). 
Key to this idea is, as discussed by Sauls, Stein & Serene (1982) (see also Vardanyan & Sedrakyan (1981)), the fact that the 
entrainment effect induces a flow in the proton fluid around each neutron vortex. This, in turn, generates a local magnetic field 
of the order of 10 14 G off of which the electrons scatter dissipatively (Alpar, Langer & Sauls 1984). The outcome is a coupling 
between the neutrons and the interpenetrating conglomerate of charged particles. This mechanism has been considered in 
two inportant scenarios. First, Alpar & Sauls (1988) have argued that it leads to the core fluids coupling on a timescale of 
400 — 10 4 rotation periods. (For alternative pictures, where the coupling timescale is significantly different, see Sedrakyan, 
Shakhabasyan & Movsisyan (1985); Sedrakian & Sedrakian (1997); Sedrakian (1998).) As this is far faster than the relaxation 
timescale following a Vela pulsar glitch, one can argue that the glitches cannot be associated with the core but rely on the 
conditions in the crust. Second, Mendell (1991b) discussed the fact that the mutual friction is also an important damping 
agent for neutron star oscillations. This is of key importance for potential gravitational-wave driven mode-instabilities (see 
Andersson (2003) for a recent review). In fact, the current thinking is that the mutual friction supresses the instability in the 
star's f-mode entirely (Lindblom & Mendell 1995). The effect on the unstable r- modes is not quite as devastating (Lindblom 
& Mendell 2000), but it could well be that the mutual friction sets the most stringent constraints on the instability window 
also in this case. 

In recent work we have applied the general formalism for superfluid neutron stars in Newtonian gravity developed by Prix 
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(2004) (see also the closely related work by Carter & Chamel (2004, 2005a,b)) to a set of problems relevant for astrophysical 
neutron stars. In particular, we have considered rotating stellar configurations where the two fluids are allowed to spin at 
different rates (Prix, Comer & Andersson 2002; Andersson & Comer 2001; Prix, Novak & Comer 2005). We have also discussed 
the nature of the inertial modes of oscillation (of which the r- modes form a sub-class) (Prix, Comer & Andersson 2004), and 
investigated the possibility that a two-stream instability may operate in a superfluid neutron star (Andersson, Comer & Prix 
2003, 2004). These papers demonstrate clearly i) that the superfluid oscillation problem is richer than tends to be assumed, and 
ii) that the entrainment effect (whereby the flow of one fluid imparts momentum in the other) plays a key role in determining 
the nature of the various modes of oscillation and the extent to which the various fluids partake in the pulsation. 

Our previous studies were based on the equations that follow after smooth-averaging over a collection of vortices. This 
"macroscopic" approach does not provide insight into the actual vortex dynamics. In order to devise a model for the required 
mutual friction force we must make connection between the large scale dynamics that we have previously considered and the 
"mesoscopic" level, which is sufficiently resolved that individual vortices can be distinguished yet sufficiently coarse that we 
do not have to worry about "microscopic" quantum effects (other than the quantisation of vorticity). Our analysis is based 
on the well-established procedure for deducing the mutual friction force in the case of superfluid Helium (Hall & Vinen 1956; 
Bekarevich & Khalatnikov 1961; Donnelly 1991; Barenghi, Donnelly & Vinen 2001), and proceeds in three main steps: First we 
discuss the nature of the quantised vorticity, making the appropriate distinction between transport velocities and momenta. 
This is important if one wants to correctly account for the entrainment. Next we derive an expression for the Magnus force, 
which describes how a bulk flow imparts a force on the vortices (analogous to the so-called Joukowski lift in standard fluid 
mechanics); see Sonin (1987) for a useful review. Finally, we derive an expression for the vortex- mediated mutual friction 
force, which couples the superfluid neutrons to the conglomerate of charged components. Having derived this expression we 
discuss the relevant coefficients and compare our final results to previous ones in the literature. 



2 THE SUPERFLUID EQUATIONS OF MOTION 

We take as our starting point the two-fluid equations derived by, for instance, Prix (2004) (see also Andersson & Comer 
(2005)). In this description the number density of each fluid obeys the continuity equation 

^ + V J (n x ^)=0. (1) 

Here, we distinguish between the "consitutent index" x which (in the present context) can be either n or p, and the spatial 
index j. The index n represents the superfluid neutrons while p corresponds to a conglomerate of all charged particles (protons 
and electrons), which are expected to flow together due to electromagnetic coupling (Mendell 1991a). In the following, repeated 
constituent indices (x and y) never imply summation while spatial indices i, j and k satisfy the Einstein summation convention. 
In Eq. (1), n x is the number density and v x is the transport velocity. That is, n x = n x « x represents the true number density 
current for species x. 

Each fluid satisfies an Euler-type equation, which ensures the conservation of total momentum. For constituent x this 
equation can be written 

{m + viv >) [v * + £xW ^ ] + Vi($ + Ax) + £xw r v< ^ = • (2) 

Here we have defined the relative velocity 

W T = v i ~ v i ■ ( 3 ) 
Furthermore, 

Mx 1 dE 

= = 7T— , (4) 

m x m x dn x 

where E is the internal energy of the system, is the relevant chemical potential per unit mass. The entrainment is included 
via the coefficients 

dE 

e x = 2p x a where a = 2 , (5) 

and $ represents the gravitational potential. For a detailed discussion of these equations, see Prix (2004); Andersson & Comer 
(2005). 

Using the fact that the momentum per particle (which is canonically conjugate to the number current n x ) of each fluid 
is defined as 

Pi = m x [vi + e x wf x ] , (6) 
and identifying the spatial derivatives in (2) as components of the Lie derivative £ Vx associated with the velocity v x , i.e. using 
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£v x Wi = viVjWi + WjVivi (7) 
which holds for a general co- vector Wi, we can rewrite (2) as 

(^ + ^ x )p" + m x Vi(* + /ix-^x) =0. (8) 

Let us now make contact with microphysics. In the case of a neutron superfluid the momentum pf is related to the 
gradient of the phase Xn of the condensate wavefunction (the "order parameter" ) via 

Pi = |v iX n , (9) 

where the factor of 2 is introduced since we are dealing with neutron Cooper pairs. Taking the curl of this relation, we see 
that the superfluid is generally irrotational. It follows from (8) that this property is conserved by the flow, see Prix (2004), 
and hence it is natural that vortices are associated with quantised momentum circulation. To some extent, this contrasts with 
the "orthodox" Landau formulation of superfluids (see the work of, for example, Mendell (1991a,b)). In that paradigm pf/m n 
is refered to as the "superfluid velocity", Vf (say). Conceptually, this is somewhat confused but there is no real risk of making 
significant mistakes as long as one does not try to account for entrainment. After all, if e n — we trivially have Vf = vf. As 
we will see, the situation when entrainment is considered is far from trivial and one must take some care in order to avoid 
inconsistencies. 

If we introduce vortices in the superfluid then it is easy to show that the circulation in the neutron fluid must be quantised. 
Integrating along a contour which encloses a single vortex we have 

C = jp*dx l = J (e^VjpDdSi = \ . (10) 

If for simplicity we assume that the vortex is straight, and introduce cylindrical coordinates (r, 6, z) centered on the vortex, 
the neutron momentum can be represented by 

Pe = l- (ID 

It should be noted that, in an orthonormal basis, this corresponds to the more familiar looking vortex solution p — Cee/2nr. 
Formally, the vorticity of the flow is associated with the singularity at r = 0. However, on a macroscopic scale we can 
meaningfully introduce a smooth averaged "rotation velocity" by comparing the standard result 

20* = e iik Vivt , (12) 

to (6). Hence, on the macroscopic level we identify (assuming that the lengthscale considered is sufficiently small that we can 
treat £ n as a constant) 

e*> k v jP l = 2m a [n* n + £n (np - K)} • (13) 

As a result, we find that by enclosing M vortices we get 

2% m n [Q n + e n (fip -OOHr = ^ , (14) 
where we have used Q l x = Q,^e\. Defining the surface vortex density n v as 

n " = 2^d~r ' (15) 
we see that 

n v H = 2[« n + £ n (a p - « n )] . (16) 

It should be noted that this result, which clearly displays the interpenetrating nature of the two fluids, differs from the 
"standard" result (see for example Alpar, Langer & Sauls (1984)). Because of the entrainment effect, the number density of 
neutron vortices depends explicitly on the rotation of the proton fluid (this point was recently discussed also by Chamel & 
Carter (2005)). 

Let us now consider the force acting on a single neutron vortex. This means that the neutron momentum is represented by 
(11), since we focus on the mesoscopic scale. Meanwhile, there should (on this scale) be no circulation in the proton momentum. 
We will motivate this assumption later when we account for the fact that the protons are charged. Taking p' vp = we have 

— =Vg + E n (v e p - Vg ) = — , (17) 
m n Z7T 



111 



= v ;p + £p (vP-vZ p ) = 0, (18) 



p 
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where k = h/2m n . For later convenience we have identified the momentum associated with the vortex flow p l m by an additional 
index v. This will later allow us to distinguish between a uniform flow past, and the rotation induced by, the vortex. If we 
solve for the individual flows induced by the vortex we find 

»" = G^tU' (i9 » 

, — • (20) 
£p + £n - 1 I 2tt V ' 



3 THE MAGNUS FORCE 

Having discussed the nature of the quantised vorticity, we are ready to investigate the force acting on a vortex due to a 
uniform flow past it. To do this, we require the flow of momentum 7r x = n x p x . That is, we need 

dnf d , „, dpf x <9n x ,„„. 

-at = «<^> = (21) 

Combining (1) and (2) we can show that 
d_ 

at 



^ (ir? + nf) = -pVi$ - n n V,/j„ - n p V;/i p + aViWp n -Vj [n n pivi + n p p?vp] 



~ -Vj {n n p>n + ripp^vp + [p$ + n n p n + n p p p - aWp n ] } 

= -V,^ , (22) 

where p = p n + p P = m n (n n + n p ) is the total mass density — we take m n = m p throughout this paper. It has been assumed 
that we are working on a sufficiently small scale that we can take the number densities n x and the entrainment parameter a 
constant. 

In order to quantify the force acting on a vortex it is natural to work in a frame in which the vortex is at rest. In that 
frame we want the background flow to be stationary and irrotational. Then we want to ask what the effects of introducing a 
single vortex in this flow may be. Translating (2) into a frame moving with velocity vf we get 

(§1 + «4V,-) [»f + e*«r] + + Ax) + e xW r V iWxL = , (23) 
where wf L — vf — vf. Imposing the condition of a stationary and irrotational flow, the latter of which leads to 
«4V 4 [w? + £xW f ] = «4Vj [w? + e x wT] , (24) 
we have 

Vi |^w xL + * + /Ux+£xW xL w7} =0 . (25) 
After integration, this yields 

iw xL + $ + /i x + £xW xL Wj y = C x = constant . (26) 
From these two equations, we have 

D = p n Cn + PpCp = ipnWnL + ^Pp^pL + + fln/Un + n p p p - 2«Wp n . (27) 



Combining this with (22), obviously translated into the vortex frame, we can show that 
III = Pnwi L w? L + p p wl >L wf-2aw pn wr + 8t 



(28) 



2 r " nL 2 f 

The force per unit length that acts on the vortex is now determined by the flow of momentum through a cylinder enclosing 
the vortex. That is, we need to evaluate 



fi = <b nisjdi , (29) 



ic 

where C encloses the vortex and Sj is the unit normal to the cylinder. For a single vortex we can assume the flow to be 
approximately of the form 

w? = U? + vr, (30) 

where Uf is uniform, stationary and irrotational in the vortex frame. The vortex flows are given by (19) and (20). 
Substituting in (28) and retaining only those terms that will eventually contribute to the integral in (29) we find 
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n& = nnUibrsj - ?T*i) + ^uiivT^ ~ p? p «0 ■ (31) 

To arrive at this result we have used the fact that all coefficients of s, which are constant on the contour vanish when integrated 
around a circle. As the vector Si in polar coordinates is proportional to the radial vector and is in the direction we 
obviously have v J m Sj = 0. Finally, as [/* is a constant flow, the integral around the vortex of {U^Sj)Uf vanishes. Our final 
expression can be rewritten as 

H-Sj = n n ti jk Ul(e klm p v i n s m ) , (32) 

where we have used the fact that p\ v = for a single neutron vortex. 

If we define the "vorticity" as a vector with magnitude k which is aligned with e ljfc Vjp£, i.e. use 

j k TTlnKi (QO\ 

£ijkPvn s — ^ ' 

(in orthonormal cylindrical coordinates k — ne z ), and work out the force integral, we arrive at the final result for the Magnus 
force acting on the vortex 

fT = P^i]kUi^ k = -p n eijkW J IlIj K k = PnCijfc^WnL • (34) 

From this result, the force density (per unit volume) on a collection of vortices follows readily as . 

It is important to note that entrainment affects the Magnus force in two ways. First of all, it enters (34) via K j . Secondly, 
it also impacts on the vortex number density n v according to (16). These results agree with the discussion of Chamel & Carter 
(2005). Our final formula (34) also agrees with the result used by Langlois, Sedrakian & Carter (1998). Mendell (1991b) makes 
use of a more generic expression which allows for the presence of vortices in each different fluid. The coupling coefficients in 
that expression are, however, left unspecified. 



4 THE MUTUAL FRICTION FORCE 

Having found the form of the Magnus force, we can determine the "mutual friction" force, which represents a balance between 
the Magnus force and standard "resistivity" due to electrons scattering off the magnetic field associated with each vortex. 
Taking the latter force to be proportional to the difference in velocity between the vortex and the charged fluid flow (protons 
and electrons), we have 

Pi = K(vf - v\) . (35) 

Assuming that the vortex can be treated as massless (see Mendell (1991a) for a justification of this assumption), this force 
must equal the Magnus force given by (34). Solving the resultant equation for v\ (using repeated cross products with k\ see 
for example Hall & Vinen (1956)) we find 



ijk 



1 



Consequently, the force per unit length acting on the vortex is 



/e = Tl{vl - vl) 
Tl 2 



P nK 2 \l+TZ 2 /pW 

,n( i 



n 2 KI + ^/pIk- 



ijk 



ijk I m lr>*7\ 

6 Kj6fc( m K W pn . (37) 



The first term in this expression is analogous to the Magnus force, although now expressed in terms of the velocity difference 
wf n . As this force is perpendicular to the relative velocity, it is non-dissipative. The second term, on the other hand, can be 
rewritten using 

e %jk K 3 e klm K l w^ = K*{tiki - g")w^ , (38) 

where the bracket can be recognised as the projection orthogonal to k l — n 1 / k. This term induces dissipation in the flow. 
From (36) and (37) we also see that: 

• In the limit TZ — > oo we must have v\ — > vf. That is, the neutron vortices are strongly coupled to the charged fluid. 

• In the opposite, "weak coupling", limit where TZ — > 0, the vortices must flow with the neutron fluid, i.e., we have v\ — vf. 

• The dissipative part of the mutual friction force, somewhat counterintuitively, vanishes in both of these limits. 
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The form of the mutual friction to be used in the equations of motion is obtained, assuming that there is no direct vortex 
interaction (see Ruderman, Zhu & Chen (1998) for a discussion of such interactions), by multiplying (37) by the vortex density 
n v . This is then the force that acts on the neutron superfluid, eg. which enters the right-hand side of Eq. (2) with x = n. An 
equal and opposite force acts on the charged conglomerate, and provides the right-hand side of Eq. (2) with x = p. 



5 ESTIMATING THE COEFFICIENTS 

To complete our investigation, and in order to facilitate the use of our results in studies of the dynamics of neutron stars, 
we need to discuss the parameters that determine the strength of the mutual friction force. In essence, we need to estimate 
the "friction coefficient" 1Z. To do this, we rely on the previous work of Alpar, Langer & Sauls (1984) and Mendell (1991b). 
Below we "translate" their analysis into our formalism. 

The contribution to the mutual friction force which is expected to provide the dominant coupling mechanism between the 
superfluid neutrons and the charged conglomerate is due to electrons scattering off the magnetic field associated with each 
vortex. An early analysis of this coupling was carried out by Sauls, Stein & Serene (1982), who discussed the importance of 
the spontaneous magnetisation of the vortex. Shortly after this analysis it was realised that the proton current induced by 
the entrainment effect would lead to a significantly stronger magnetic field (Alpar, Langer & Sauls 1984). Hence, we focus 
our attention on this case. 

For superconducting protons, the expression for the momentum must be replaced by (see, for example, Prix (2005)) 

V\ + \M = ^V iXp • (39) 

(Here and in the following we are using Gaussian units). The magnetic field associated with the flow follows from, firstly the 
definition of the magnetic potential 

Bi = e rjk V 3 A k (40) 
and secondly, the Maxwell equation 

ji = ^e ijk V j B k = en p vf , (41) 

where the right-hand side is the charge current. In writing down this relation we have adopted the convention that the charge 
currents affect only the magnetic induction Bi, see Tilley & Tilley (1990) for further discussion. It is now straightforward to 
combine (40) and (41) to obtain the London equation for Ai. 

Presently, we are primarily interested in the magnetic field generated by the flow of entrained protons around a single 
neutron vortex. This means that it is natural to assume that there are no fluxtubes in the proton fluid. This is equivalent to 
assuming that the phase x P is smooth. Then a gauge transformation can be made such that (39) is replaced by, see Tilley & 
Tilley (1990), 

Pi + -^Ai = . (42) 
Thus it follows that 

m p (l — £ P )vf = —-Ai — m p e p Vi , (43) 

and we arrive at the following equation for the magnetic flux (using the fact that X7iB l — and assuming that the vortex is 
represented by a delta-function) 

Here, the effective London penetration length A, is defined by 
2 _ cjm| 2ap\ ( \ 2a 



M = 1 - ^ 1 - — . (45) 

4ne 2 p p \ p n p p J \ p n J 

Comparing (44) to equation (14) of Alpar, Langer & Sauls (1984), and recalling the fact that the "superfluid velocity" in the 
orthodox formalism is in fact the rescaled momentum, we find that we should identify 

^ P = p P fl-^y7l--) , and pT = -2a(l-^Y . (46) 

V PnPp J \ pn J \ PnPp J 

The relationship between the two formalisms has already been discuseed by Prix, Comer & Andersson (2002). Using their 
analysis we readily demonstrate that the above identification is correct. 

As the equations for the magnetic field are identical, the required solution is identical to that given in previous work 
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(Alpar, Langer & Sauls 1984; Mendell 1991a). Solving (44) for the case where the vortex is represented by a delta-function ( 
Fetter & Hohenberg 1969), we readily find that the only non- vanishing component of the magnetic field is 



B 



rK (r/A t 



27rA2 J 

Here Ko is a modified Bessel function, 

he m p pv n 
2e m n p p s 



_ iil nip fj s 

^* n_ -PP ' 



and m p = m n for all practical purposes. Crucially, this means that 



for r » A» 



\A\ cx Ki{r/K,) ~ \ — exp(-r/A*) 



(47) 



(48) 



(49) 



This resolves the apparent contradiction between (18), which formed a key part of our derivation of the Magnus force, and the 
correct formula for charged protons, Eq. (42). We see that the analysis in Section III holds, provided that the force calculation 
can be performed sufficiently far away from the vortex that Ai can be neglected yet close enough that the assumption of 
essentially constant densities, entrainment parameters etcetera holds. This should always be possible, since A* is many orders 
of magnitude smaller than eg. the intervortex separation. 

For later convenience, it is useful to take a brief detour at this point and introduce the effective proton mass m*. This 
concept was also discussed by Prix, Comer & Andersson (2002) . It is natural that the entrainment effect can be expressed in 
terms of an altered "effective" mass since it couples the momenta of the two fluids. The analysis leading to this is, in fact, 
very simple. In a frame comoving with the neutrons, i.e. in which vf = 0, (and well away from all vortices in a region where 
we can neglect the magnetic field, according to the discussion above) we have 



* p 
m p v i 



m p (l - e p ) 
Hence, it follows that 

2a = p p £ p = n p (m p - m*) = n p <5m* 



(50) 



(51) 



(Note that we could alternatively have defined the effective mass in the frame where pf = 0. As discussed by Prix, Comer & 
Andersson (2002), the result is the same provided that the proton fraction is small — the limit which is relevant for neutron 
star cores.). 

This means that we can use 

-l 



np 

Ps = -Pp 



Snip 



1 - 



5m* 



-Pp 



Snip 



which should be accurate in the case of neutron stars where the proton fraction x p — p p /p is small. We also have 



pp P n m p - Ppdnip 
Ps = P P - 



p n mp — pSmp 



: Pp 



(52) 



(53) 



These expressions are, not surprisingly, identical to those given by Alpar, Langer & Sauls (1984). Given these approximations, 
the penetration length is 

1-1/2 



A, « 1.3 x 10 



1 0.05/ 



10 14 g/cm 3 



fm 



and the magnetic field associated with the vortex core is approximated by (Alpar, Langer & Sauls 1984) 



B : 



1**1 
2ttA2 



1.9 x 10 14 G 



1 0.05/ 



lO^g/W 



8m* 



(54) 



(55) 



In order to estimate the relaxation time for electrons scattered off the vortex magnetic fields, we combine three further 
results from Alpar, Langer & Sauls (1984). The first is the relaxation timescale to in the limit of a vanishing vortex radius. 
It follows as 



t 1 = ivN T $l , 
where 

N - —n ( eH V 

% v \2m c c/ 



Efc 



Efh I (tyHc) 2 



(56) 



(57) 



We will assume that the electrons are ultrarelativistic, i.e. use 
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E Fc = hck Fc , where k Fc = (37r 2 n ) 1/3 = (37r 2 n p ) 1/3 . (58) 
As the electrons and protons are expected to couple on a much shorter timescale we account for the increased inertia by using 



_ / m p c 2 \ 
I hck Fc J 



T() . 



(59) 



The final factor encodes the dependence on the finite size of the scattering centre. As discussed by Alpar, Langer & Sauls 
(1984), this leads to 



16 a 



where 



a 
1 



2fc Fo A» ?a 120 



1/2 r 



( x p 
V0.05 



Pp 



10 14 g/cm ; 



-1/6 



(60) 



(61) 



Having arrived at an estimate of the timescale on which the vortices relax to the motion of the charged components, we 
can make connection with our expression for the mutual friction force from the previous section. To do this we note that the 
relative velocity between vortices and charged components relaxes according to dtAvi = —Avi/r v , from which we can deduce 
that the average force (per unit length) acting on a typical vortex is 



(fi) = 



Pp 



7lv Tv 

Comparing this to (35) we see that 



Pp 



(62) 



(63) 



« 1.6 x 10" 7 ] 
pn/t I \ m v \mt V0.05/ 



As discussed in the previous section, it is useful to establish whether we are in the regime of strong or weak coupling. The 
above analysis leads to the estimate 

,c *X4 , x 7/3 ( v 1/3 

% « 1 ■ (64) 

v 10 14 g/cmV 

Given that the effective proton mass is such that mj/m p w 0.5 — 0.7 we are firmly in the weak coupling regime. 

Finally, let us return to (37) — the force per unit length acting on an individual vortex — and replace it with the 
relevant force acting on the superfluid neutrons after averaging over the vortices. This is a slightly subtle issue. First we need 
to appreciate that the vortices are already accounted for in the averaged equations, eg. (8). In a sense this means that the 
Magnus force (34) is already contained in this equation. Hence we need to add only the resistive part (35) to the description. 
The force that we require thus follows simply by multiplying (37) by the local surface density of vortices n v . Provided that 
we are dealing with the weak coupling limit, the resultant force acting on the neutron fluid can be written 

/i k> iik ~ I m . k>' ilk pn 

m { = Bp n n v e J KjCkimK w pn + B p n n v e J K-jW^ . (65) 



An equal and opposite force acts on the proton fluid. Here 

S=A,4x 10- ( S <] 2 (^.) (^) 7/6 ( 
p n n y m p J ymp J V0.05/ y 

and 

B' = B 2 . 



10 14 g/cm d 



1/6 



(66) 



(67) 



are dimensionless parameters. These results should be compared to the parameters used by Mendell (1991b), and it is easy 
to confirm that the two results are in perfect agreement. 



6 CONCLUDING REMARKS 

In this paper we have derived the form of the vortex-mediated mutual friction, which arises as electrons scatter dissipatively off 
of the magnetic fields associated with the entrained proton currents and each neutron vortex, within the superfluid formalism 
developed by, for example, Prix (2004); Andersson & Comer (2005). In doing this we have made contact with previous 
work based on the orthodox Landau formulation (Alpar, Langer & Sauls 1984; Mendell 1991a), and demonstrated that the 
two pictures are consistent. Since our description incorporates the entrainment effect in a transparent way (by making the 
appropriate distinction between true transport velocities and momenta) this comparison lends strong support not only to our 
present results but to the previous work as well. 
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To conclude our discussion, let us put the final expression for the mutual friction force (65) to use by working out the 
timescale on which a difference in rotation between the two fluids in a neutron star core (eg. following a pulsar glitch) is 
relaxed locally. Given (8) and (65) we see that the system evolves according to 

n n dtPi + ■ ■ ■ = /f f 1 m* np BkUv 



ripOtpV + ... = -/, 



mf 



^ f < p + ...«-=< p , (68) 



where, on the left-hand side we have used the definition of the momenta and assumed that the proton fraction is small 
(in accordance with the preceding analysis), and on the right-hand side we implicity assume that the velocity difference is 
perpendicular to the vortices (that is, the two fluids rotate around the same axis). From this expression we see that the 
timescale on which the two fluids are dynamically coupled can be estimated by 

T d *^jp^. (69) 
m p Ben, 

Taking n v K « 4-k/P, i.e. assuming that the two rotation rates are similar to the observed pulsar period P, cf. (16), we have 

-1/6 / \ - 1 / 6 



*" 10P «U*J {^) U^wJ • (70) 

This estimate is about one order of magnitude smaller than the classic result of Alpar & Sauls (1988). That the two results 
differ is perhaps not too surprising. After all, the Alpar & Sauls (1988) analysis was based on an explicit solution for the 
motion of an individual vortex affected by the Magnus force (34) and the resistivity (35) given a constant relative rotation 
rate. In contrast, our estimate does not refer to the explicit vortex motion, only to the way that the two fluids couple via the 
vortices. Of course, the main astrophysical conclusion remains unaltered. The coupling timescale is much shorter than the 
observed relaxation timescale following (say) the large Vela glitches (Alpar & Sauls 1988). This suggests that the glitches are 
unlikely to be associated with the core and points instead to the superfluid neutrons in the crust playing a key role. 

Basically, we have now prepared the ground for discussions of the relevance of mutual friction in different astrophysical 
scenarios within our formalism. This is a very important step forwards since it allows us to consider key problems concerning, 
for example, the mutual friction damping of pulsation modes driven unstable by gravitational radiation. As the answer may 
be of significance for gravitational-wave observations, the available results in that problem area (Lindblom & Mendell 1995, 
2000) must be verified by independent work. Given the present analysis, we are set to carry out such calculations and expect 
to report the results in the not too distant future. 
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